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A possible solution for the problem of memory-size and computer-time, is 
the extrapolation of basis-set^. This extrapolation has two exponents a and 
(5, corresponding to the HF (reference energy) and the energy of correlations 
(EC), respectively. For a given system, the exponents are taken as constant^, 
and potential energy surfaces (PES) are generated. We have found that the 
values of a and (5 are not constant, but vary from position to position in 
the system. How to deal with such situation and get very accurate PES, is 
discussed. 

Truhlar^ proposed a nice concept for extrapolation of basis-set, where the total energy 
is expressed as the sum of HF and EC. In the method, two exponents a and I3, respectively, 
are introduced. 

The calculations based on the basis-sets are required in all branches of science and 
technology. In the Astronomy &: Astrophysics, their contribution is remarkable. It plays 
important role for analyzing spectrum, coming from the interstellar medium (ISM). 

For getting information about far distant cosmic objects, one depend on the radiations 
lying in the, which are generated by the molecules, present in the object. The kinetic 
temperature in an object having molecules is few tens of Kelvin, in general. So, in most 
of the cases, scientists are concerned with the pure rotational levels. 

We consider the molecule thioformaldehyde (H 2 CS), for example, whose identification^ 
in the ISM was considered as an achievement, because earlier attempts made several times 
remained unsuccessful, though they were based on very accurate laboratory studies. For 
analyzing spectrum, one considers an appropriate number of energy levels. These levels 
are connected through radiative and collisional transitions. Galculation of collisional cross 
sections is the most difficult task in the investigation. But, the scientists are interested to 
have accurate cross sections. 

The collisions are due to the most abundant molecule H 2 . However, many scientists 
do not consider internal structure of H 2 , and therefore they replace H 2 by the He atom, 
as both the H 2 and He have two electrons and two protons, and the interaction potential 
depends on the charges of the interacting particles. For the calculation of collisional 
cross sections, the interaction potential between H 2 CS and He is required, for which the 
scientists are putting their best efforts. 

The interaction potential is calculated with the help of MOLPRO'^ by using the best 
method CCSD(T)^ and the best basis-set the aug-ccpVnZ, where n denotes the number 
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of zeta used and it stands for one of the D (Double-zeta), T (Triple-zeta), Q (Quadruple- 
zeta), 5 (Quintuple-zeta), etc. The prefix ’aug’ denotes the augmented versions of the 
basis set. With the increase of n, the accuracy of results increases, but the requirement of 
memory-size and computer-time increases exponentially, which is not affordable, beyond 
n = 6 (say). 

How to reach n = oo by using the values for n = 2 and n = 3 is a nice idea used 
in the basis-set extrapolation^. That is, one has to calculate the HF and EC by using 
the CCSD(T) method along with the aug-ccpVDZ and aug-ccpVTZ basis-sets, separately. 
With the help of the energies of these two sets, the limiting energies (corresponding to 
n = oo) can be obtained. For a given system, the values of the exponents a and (3 were 
taken constant. 

We put a question. Are the values of a and (5 constant through out the system? The 
answer to this question can be found in the following manner. We prepare Figure 1, where 
we plot F[= (4“'>' — 2 ~^)/[ 2 >~^ — 2“''')] versus 7 . It is a nice smooth curve. 

The position of He atom relative to the center-of-mass of H 2 CS is expressed in terms 
of the spherical polar coordinates (r, 9, 4>). The interaction potential between the H 2 CS 
and He is calculated for various positions (r, 6, ip). 

We calculate HF and EC for one more basis-set, aug-ccpVQZ, i.e., n = 4. Let the values 
obtained by using the basis-sets aug-ccpVDZ, aug-ccpVTZ, aug-ccpVQZ, be denoted by 
1 ^ 2 , V 3 , V 4 , respectively. We calculate the parameter F = (V 4 — V 2 )/(h 3 — V 2 ), and from 
Figure 1, find out 7 , corresponding to this value of F. The 7 is either a or /3, corresponding 
to the HF or EC, respectively. For better accuracy, this graphical method is replaced by 
the numerical interpolation. In this way, we find the values of a and /3 at each position. 

We have considered 532 positions specified by (r, 6, (p), with r = 2 (0.25) 6.5 A, 9 = 
0° (30°) 180°, (p = 0° (30°) 90°. The obtained values of a and /3, are plotted in Figures 
2 and 3, respectively, as a function of position. Figures 2 and 3 show that neither a nor 
13 is constant. Both of them vary from position to position in the system. For very large 
r, the a and /3 attain the limiting values, denoted by ao and /Iq. The values obtained are 
ao = 3.157 and /3o = 2.548. These are the values which are obtained with the help of 
statistical methods and are in practice to use. For small r, the deviations of a and f3 with 
respect to their limiting values, ao = 3.157 and f3o = 2.548, are very significant. With the 
increase of r, the deviations decrease. 

In the basis-set extrapolation, the limiting value 14o {n = 00 ) is expressed as Foo = 
AV 3 — {A — 1 )^ 2 , where A = /{3F — 2'^). Using the constant values ao = 3.157 and 

/3o = 2.548, we have calculated total energy Eo at each position. For the values of a and 
p3 at each position, we have calculated total energy E, and plotted in Figure 4, the value 
of {E — Eo) as a function of position. The values of {E — Eq) are up to the order of 10“^ 
atomic unit (A U). The deviations shown in Figure 4 are quite significant, in particular 
for small r, as the total energy is calculated up to the accuracy of 10“® A U. 

Finally, by using these values of a and /3 at each position, the HF and EC, and finally 
the total energy, corresponding to n = 00 , can be calculated, and very accurate PES can 
be generated. 
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Eigure 1; Variation of F [= (4“'’' — 2“^)/(3“'>' — 2“'*')] versus 7 . The F is plotted on the horizontal 
axis and 7 on the vertical axis. 


Eigure 2; Values of a for various positions. We have r = 2 (0.25) 6.5 A. Each r has 28 positions. 


Eigure 3: Values of j3 for various positions. We have r = 2 (0.25) 6.5 A. Each r has 28 positions. 


3 



Figure 4: Values of {E — Eq) for various positions. We have r = 2 (0.25) 6.5 A. Each r has 28 
positions. 
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